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0-^. Abstract 

< 

Sand pile models are dynamical systems describing the evolution from N stacked grains to a 

stable configuration. It uses local rules to depict grain moves and iterate it until reaching a fixed 

configuration from which no rule can be applied. Physicists L. Kadanoff et al inspire KSPM, 

extending the well known Sand Pile Model (SPM). In KSPM(D), we start from a pile of A^ stacked 

grains and apply the rule: D—1 grains can fall from column i onto columns i + 1, i + 2, . . . , i + D—1 

if the difference of height between columns i and i+1 is greater or equal to D. Toward the study of 

Y^ \ fixed points (stable configurations on which no grain can move) obtained from A^ stacked grains, we 

propose an iterative study of KSPM evolution consisting in the repeated addition of one grain on 
a heap of sand, triggering an avalanche at each iteration. We develop a formal background for the 

^ \ study of avalanches, resumed in a finite state word transducer, and explain how this transducer 

^\ ' may be used to predict the form of fixed points. Further precise developments provide a plain 

formula for fixed points of KSPM(3), showing the emergence of a wavy shape. 



Keywords: Discrete Dynamical System, Self-Organized Criticality, Sand Pile Model, Fixed 
point, Transducer. 



1. Introduction 



X 

H ■ This paper is about cubic sand grains moving around on nicely packed columns in one dimension 

- - - (the physical sand pile is two dimensional, but the support of sand columns is one dimensional). 

We follow the arbitrary convention that when sand grains can only fall in one direction according 

to iteration rules, this direction is the right. So when there is no ambiguity, we use any variation 

of the words right and left to refer to the direction of grain falls and its opposite. 

1.1. The framework 

Sand pile models were introduced in [l| as systems presenting a critical self-organized behavior, 
a property of dynamical systems having critical points as attractors. In the scope of sand piles. 
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starting from an initial configuration of N stacked grains the local evolution of particles is described 
by one or more iteration rules. Successive applications of such rules alter the configuration until 
it reaches an attractor, namely a stable state from which no rule can be applied. Self- organized 
criticality (SOC) property means those attractors are critical in the sense that a small perturbation 
— adding some sand grains — involves an arbitrary deep reorganization of the system. Sand pile 
models were well studied in recent years. 

Starting from N stacked grains on column and no grain elsewhere, the first one dimensional 
Sand Pile Model (SPM) applies the rule: if the difference of height between columns i and z + 1 is 
greater or equal to 2, then one grain falls from column i to column i -\- 1. In [9| the authors show 
that the set of reachable configurations endowed with the successor relation given by the iteration 
rule is a lattice. They furthermore provide a simple characterization of reachable configurations 
and give a plain formula describing the unique fixed point according to the number N of grains. 
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Finally they prove the convergence time to be in 0{n2). For a survey on SPM, see 10(]. A slight 
variant of SPM has been studied in [J]: instead of applying the rule sequentially — once at each 
iteration step — , the rule is applied in parallel on all possible positions. This model is called Parallel 
SPM (PSPM) and is deterministic. It reaches the same fixed point as SPM, but in time 0{n) (very 
smart and technical proof in [J]). [5(] and [20(| explore Symmetric SPM (SSPM) where grains can 
fall either to the right or to the left according to symmetric rules. The authors both give a simple 
characterization of fixed points shapes (there is no lattice structure anymore) in this sequential 
model. The Parallel SSPM (PSSPM)jWhich is non deterministic since choices may occur on the 
top column, has been considered in p where a constructive proof shows that the fixed points 



shapes are the same as in SSPM, and in 16|] where the authors compare not only shapes, but also 
positions. 

Another interesting question about SPM is the prediction problem (namely, the problem of 
computing the fixed point), which has been proved in [l5| to be in NC'^ for the one dimensional case 
(it means that the time needed to compute the fixed point is in 0(log N) on a parallel computer 
with a polynomial number of processors, where N is the number of grains), and P-complete when 
the dimension is > 3. 

In [131) Kadanoff proposed a generalization of classical models closer to physical behavior of 
sand piles, in which more than one grain can fall from a column during one iteration. Informally, 
Kadanoff sand pile model with parameter D, KSPM(D), is a discrete dynamical system whose 
initial configuration is composed of a finite number N of stacked grains, moving in discrete space 
and time according to an iteration rule: if the height difference between column i and i + 1 is 
greater or equal to D, then D — 1 grains can fall from column i to the D — 1 adjacent columns on 
the right (see figured]). Note that KSPM(2)=SPM. 



> D 




Figure 1: KSPM(Z)) iteration rule. 
Figure [2] presents an example of evolution. 
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Figure 2: Example of an evolution from 24 stacked grains to the associated stable configuration. 
At each step, the arrow is labelled by the index of the column on which the rule has been applied. 



1.2. Definitions and previous results 

More formally, sand pile models we consider are defined on the space of ultimately null de- 
creasing integer sequences. Each integer represents a column of stacked sand grains and iteration 
rules describe how grains can move from columns. Let h = {ho, /ii, /12, • • • ) denote a configuration 
of the model, where each integer hi is the number of grains on column i. Configurations can also 
be given as sequences of height differences a = (ctq, cji, cj2, . . . ), where for all i > 0, ai = hi — /ij+i. 
We will use this latter representation throughout the paper, within the space of ultimately null 
non- negative integer sequences. 

Definition 1. The Kadanoff sand pile model with parameter D, KSPM(D), is defined by: 

• A set of configurations, consisting in ultimately null non-negative integer sequences. 

• A set of iteration rules.- we have a transition from a configuration a to a configuration a' on 
column i, and we note a ^ a' when 



- cr. 

- cr. 

- cr; 



/ 



= C7i_i +D-1 (fori^O) 
(Ji-D 



'i+D~l - (^i+D-l + 1 
- cr'j = CTj for j ^ {i-l,i,i + D - 1}. 

We also say that i is fired. 

Remark that according to the definition of the iteration rules, a condition for a' to be a con- 
figuration is that ai > D. We note a ^- a' when there exists an integer i such that a ^ a' . The 
transitive closure of — )■ is denoted by — )■, and we say that a' is reachable from a when a ^ a' . 

A basic property of the KSPM model is the diamond property. If there exists two distinct 

integers i and j such that a ^ a' and a — t- a" , then there exists a configuration a'" such that 
a' 4 a'" and a" A a'". 
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Figure 3: The set of reachable configurations for A^ = 24 and Z) = 3. The initial configuration is 
on the top, and for short we denote 0^ the infinite sequences of 0. 



We say that a configuration a is stable, or a fixed point if no transition is possible from a. As a 
consequence of the diamond property and the termination of the evolution on finite configurations, 
one can easily check that, for each configuration a, there exists a unique stable configuration, 
denoted by 7r(cr), such that a — )• Tr{a). Moreover, for any configuration a' such that a — )• a', we 
have 7r{a') = ir^a) (see [I4] for details). 

We are interested in the evolution from a finite number N of grains to a stable configuration. 
Figure [3] depicts the set of reachable configurations for D = 3 from the initial configuration with 
24 grains. Applying the rule once on the initial column leads for example to (21, 0, 1, 0, 0, 0, 0, . . . ) 
and we denote 0"^ the infinite sequence of 0, hence we write (21,0,1,0"^) that configuration. For 
short, we state 7r(iV) = 7r((A^,0'^)). 

Sand pile models are specializations of Chip Firing Games (CFG). A CFG is played on a 
directed graph in which each vertex v has a load l{v) and a threshold 9{v) = deg~^{v]LJ, and the 
iteration rule is: if l{v) > 6{v) then v gives one unit to each of its neighbors (we say v is fired). As 
a consequence, we inherit all properties of CFGs. 



(v) denotes the out-degree of v. 



Kadanoff sand pile is referred to as a linear chip firing game in 1^. The authors show that 
the set of reachable configurations endowed with the order induced by the successor relation has 
a lattice structure, in particular it has a unique fixed point. Since the model is non-deterministic, 
they also prove strong convergence i.e., the number of iterations to reach the fixed point is the 
same whatever the evolution strategy is. The morphism from KSPM(3) to CFG is depicted on 
figure m 

When reasoning and writing formal developments about KSPM, it is much more convenient 
to think about its CFG representation because it is independent of the height of columns. Let us 
recall that throughout the paper, we consider sequences of height differences (except when explicitly 
specified) and the associated iteration rules where units of height difference move between columns. 




Figure 4: The initial configuration of KSPM(3) is presented as a CFG where each vertex corre- 
sponds to a column (except the sink, vertices from left to right corresponds to columns 0, 1, 2, 3, ... ) 
with a load equal to the difference of height between column i and i + 1. For example the vertex 
with load N is the difference of height between column {N grains) and column 1 (0 grain). 



A recent study ([lll|) showed that in the two dimensional case the avalanche problem (given a 
configuration a and a column i on which we add one grain, does it have an influence on index j?) 
on KSPM(L') is P-complete, which points out an inherently sequential behavior. 

1.3. Our contribution 

The aim of this paper is to study the fixed points of the Kadanoff Sand Pile Model, and to 
describe them with a plain formula according to the parameter D and the number of grains N. In 
section [2] we propose an inductive construction of the fixed point for a fixed parameter D and N 
grains according to the fixed point with A^ — 1 grains: we add an N^^ grain on column of the 
fixed point with A^ — 1 grains, which triggers an avalanche, ending when the fixed point with A^ 
grains has been reached. An example of avalanche is given on figure El We study the process of 
avalanches, and its emergent regularities. The main result of this section is a precise description of 
avalanches (Theorem [T|) when a property called density is fulfilled. We also prove for D = 3 that 
on the right of a column n in 0(log A^), the A^ first avalanches have this property (Proposition [2]), 
therefore we understand the behavior of the A^ first avalanches for D = 3, starting from column n. 

Prom those regularities, we build words called traces, recording the way any grain has ever 
crossed a particular column n during the A^ first avalanches. Section [3] outlines the construction of 
a finite state word transducer (Lemma HD, outputting the trace on column n + {D — 1) from the 
trace on column n. i iterations of this transducer allows us to predict how any grain has crossed 
column n + i{D — 1), hence the behavior of avalanches from column n + i[D — 1) (see figure [5a|) . 
From the trace on a column n + i[D — 1) for the A^ first avalanches, we can straightforwardly 
compute the fixed point reached from A^ stacked grains on the right of column n + i{D — 1). For 
Z? = 3, we also prove that for a number of iterations i in 0(log A^), the trace on column n + i{D — 1) 
is a prefix of (ab)'^ (Corollary [1]). Consequently, regularities on traces induced by the transducer 
let us study the asymptotic form of fixed points on the right of column n + i{D — 1). 
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At the end of each section, we give a detailed study of the case D = 3, ending in an asymptotic 
description of its fixed points according to the number N of grains, pictured on figure [6bl Unfor- 
tunately, our complete study for the case D = 3 does not trivially generalizes to any parameter D. 
Nevertheless, we hope that the present work already conveys interesting ideas on the handling of 
complex behaviors in discrete dynamical systems. 



Some partial preliminary versions of the present work previously appeared in IT]] and 18l ]. 
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Figure 5: An example of avalanche: starting from vr(24), we add one grain on column (darkened 
on the leftmost configuration) and apply the iteration rule until reaching 7r(25). At each step, the 
arrow is labelled by the index of the fired column. 

Describing and proving regularity properties, for models issued from basic dynamics is a present 
challenge for physicists, mathematicians, and computer scientists. There exists a lot of conjectures, 
issued from simulations, on discrete dynamical systems with simple local rules (sandpile model [31 
or chip firing games, but also rotor router |l4i] . the famous Langton's ant [3] [a]...) but very few 
results have actually been proved. 



2. Avalanches 

In order to study the fixed point for a parameter D and N grains, we will study the sequence 
of fixed points for a parameter D and no grain, then D and 1 grain, D and 2 grains, D and 3 
grains, etc... Throughout the paper, the parameter D is supposed to be fixed and can take any 
value greater or equal to 2, except when a particular value is specified. Computing the fixed point 
with k grains, given the fixed point with k — 1 grains, is very simple: we add one grain on column 
and perform all the possible transitions. We call this process an avalanche. 

We first explain why adding a grain on column of the fixed point with k — 1 leads to the 
fixed point with k grains in subsection 12.11 and then formally define an avalanche. If an avalanche 
verifies a certain property (if it is dense), then we show that its evolution is somehow linear 
(subsection 12. 2p , and we give a precise description of its mechanism by the mean of particularly 
unstable columns named peaks (subsection 12. 3p . The linear process description holds only when 
avalanches are dense, which happens intuitively and experimentally starting from a very small 
index, but we only managed to prove for D = 3 that avalanches are asymptotically completely 
described as linear processes (subsection 12. 4p . Asymptotically completely means that if we consider 
the avalanche occurring from the fixed point with A*" — 1 grains to the fixed point with N grains 
for a parameter D, denoting size{D, N — 1) the size of 7r(A^ — 1) (number of non empty columns) 
and C{D,N) the column index starting from which the avalanche is dense (the property is then 



true from column C{D, N) to column size{D, N — 1) 
N tends to +oo. 



1), then the ratio 



C(D,N) 
size{D,N-l) 



tends to as 



irregular i regular 
avalanches ' avalanches 




irregular 
avalanches 



reeular 
avalanches 



? e(vW) 

(a) Avalanches become regular very quickly. 
The regularity allows us to capture the way 
grains behave in a word u. u describes com- 
pletely how grains cross a particular column 
(the long vertical dashed line on the figure). 
u barely consists in the concatenation of the 
relative positions of unstable columns, called 
peaks, for the A'^ first avalanches, within the 
D — 1 columns preceding the vertical dashed 
line. Avalanches regularity also allows us, 
knowing the word u for a column i (the dashed 
vertical line), to compute the word v describ- 
ing how grains cross the column i + D — 1 
(where D is the parameter of the model). The 
computation of v is made from u via a fi- 
nite state word transducer tC, which n*'' itera- 
tion computes the word describing how grains 
cross column i + n{D — 1). 




O(logAf) 
O(logAf) 

(b) From the knowledge of the word u de- 
scribing how grains cross a column i, it is 
very easy to compute the configuration for ev- 
ery column j > i (the concatenation of every 
words we obtain on a particular column de- 
scribes how any single grain has ever cross 
that frontier). A precise study of the trans- 
ducer for D = 3 outlines, starting from any 
word, the very quick emergence of periodic 
words describing how grains cross columns 
as we iterate the transducer C 0(log N) it- 
erations of ^ outputs periodic words of the 
form ababab .... Finally, The regularity of 
the words involves a regularity on the fixed 
points, which have asymptotically a com- 
pletely wavy shape 212121 .... 



Figure 6: Presentation of the method ()6ap and its application to KSPM(3) (J6b 



2.1. Inductive construction 

In this paper, we are interested in the iterative process defined below. Starting with no grain, 
we successively add a single grain on column 0, and make all the possible firings until a fixed point 
is reached. 

Let o" be a configuration, a^^ denotes the configuration obtained by adding one grain on column 
0. In other words, if o" = (o"o,cTi, ...), then o"-'''' = (cro+l,cTi, ... ). Let A; > 0. Remark that 7r(A;—l)~'''' 
is a reachable configuration from (fc, 0'^), because we can use any sequence of firing from [k — 1, 0"^) 
to 7r(A; — 1) to get an evolution from (A;,0'^) to 7r(A; — l)-'-'^. Thus, with the uniqueness of the fixed 
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point reachable from {k,0^), we have the recurrence formula (see figure [7]): 

TT{Tr{k - 1)^°) = Tr{k) 
which, with the initial condition: 7r(0) = 0*^ allows an inductive computation of iT{k). 



I 



^(73) 



7r(73)^o 



^(74) 



Figure 7: D = 3, adding a grain on column of vr(73) gives 7r(73)~'''' and performing all the possible 
transitions until reaching a fixed point leads to 7r(74). 



2.2. Avalanches and peaks 

We are now interested in the description of the evolution from 7r{k — 1)^^ to 7r{k). However, 
from the non-determinacy of our model, this evolution is not unique. To overcome this issue we 
give below a formal definition of an evolution (a strategy) and distinguish a particular one from 
7r(A; — 1 )''''' to 7^{k), which we think is the simplest, and define it as the k avalanche. 

A strategy is a sequence s = (si, . . . , st)- We say that a' is reached from a via s when a -4 
a" -^ . . . -^ a' and we note o" — )• o"'. We also say, for each integer t such that < t < T, that the 
column St is fired at time t in s (informally, the index of the sequence is interpreted as time) . 



i}. Let s 

— i„i 







s be two 



<^cr^ 



a 



For any strategy s and any nonnegative integer i, we state \s\i = #{t\st - 

s n s 1 r I p 

strategies such that a ^ a and a ^ a . We have the equivalence: [ V i, |s 
A strategy s such that cr —t- cr' is called leftmost if it is the minimal strategy from a to a' according 
to lexicographic order. A leftmost strategy is such that at each iteration, the leftmost possible 
transition is performed. 

The k avalanche s is the leftmost strategy from 7:{k — l)^^ to 7r{k). Toward the study of fixed 
point configurations, we first consider the process of avalanches. Informally, we want to describe 
what happens when a new grain is added in a previously stabilized sand pile. For D = 2, i.e., the 
classical SPM, this description is easy: the added grain moves rightwards until it reaches a stable 
position on two consecutive columns of same height. But, for D > 2, the situation is not so simple. 

The first Lemma is rudimentary but allows to simplify some notations for the rest of the paper. 



Proposition 1. 

|s|iG{0,l}. 



For each strategy s such that it{N)^^ — )• 7r(A^ + 1) and any column i £ N, we have 



Proof. Let s = (si, . . . , st) be a strategy such that 7r(A^)~'''^ — ;■ 7r(A^ + 1). We have to prove that, 
for 1 < / < 771 < r, we have s/ ^ Sm (obviously, \s\i > for all i). To do it, we prove by induction 
on t < T that for 1 < I < m < t, we have si ^ Sm- 

For initialization this is obviously true for t = 1. Now assume that the condition is satisfied 
for an integer t such that t < T, and let z be a column such that there exists an integer / < t such 
that i = si. Let a be the configuration such that 7r(A^)~''*' -A- ... -^- a. 

Notice that the transitions which can possibly change the value of the current configuration at 
i could be: i (which decreases the value by D units), i + 1 (which increases the value hy D — 1 
units) or i — D + 1 (which increases the value by 1 unit). 

Thus we have Ui < vr(A^)| — D + D — l + l since by definition, between 7r(A^) and a, exactly 
one transition has occurred in i, at most one transition has occurred in z + 1, and at most one 
transition has occurred in i — D + 1. For i > 1, we get Uj < Tr{N)i. On the other hand, since 7r(A^) is 
a fixed point, we have: iT{N)i < D, which guarantees that st+i ^ i. For i = 0, there is no possible 
transition in i — D + 1, thus we get o"o < 'k{N)q — D + D — I, which is ctq < 7r(A^)o + 1 — D + D — I. 
Thus (To < tt{N)q < D which also gives: st+i ^ 0. 

This ensures that the result is true for t + 1, and, by induction, for T. EH 

When talking about an avalanche s. Proposition [Hallows us to write i G s and i ^ s instead of 
\s\i = 1 and \s\i = since no other value is possible. We denote by S[j j] the subsequence of s from 
i to j included. 

We first explain the "pseudo locality" of avalanches: at a time t + 1, a fired column can't be at 
distance greater — neither on the left nor on the right — than D — 1 of the greatest fired column 
of S[i,t]- Imagine, during an avalanche, that you follow the greatest fired column with a frame of 
size 2{D — 1) — 1 centered on it. Lemma [1] tells that you won't miss any firing. 

Lemma 1. Let s^ = (sj^, . . . , s^) be the k^^ avalanche. Let rt = maxs^^^ ^, . 

• Assume that s^_^_i < rt. Then s^^-^ is the largest column number satisfying this inequality, 
which has not yet been fired at time t. In other words: 

s^j^i = max{i \ i <rt and i ^ Sm ^i} 

Moreover, we have rt — s^^i < D — 1. 

• Assume that St^^ > rt. Then we have s^_,_]^ < rt + D — 1. 

Proof. We order fired columns by causality. Precisely, a column i has two potential predecessors, 
which are i + 1 and i — D + l. State i = s^. These columns are predecessors of i if they are elements 
of Sm 1, i.e if they are fired before i. Using the transitive closure, we define a partial order relation 
(denoted <caus.) on fired columns for s^. 

Now, consider the set Aj+i of ancestors of St^i {i.e. the set of columns i such that i <caus. Sf+i) 
and the set St of columns which have s^ as an ancestor {i.e. columns i such that s^ <caus. i)- We 
necessarily have rt E ^t+i. Otherwise, we have At+i Pi S'i = 0, and this allows another strategy s', 
constructed from s^ by postponing the transitions at rt and elements of St after the transition on 
Stj^i. This contradicts the fact that s is leftmost. 

Let (io, ii, . . . , ip) be a finite sequence such that io = rt, ip = St^i and, for each j with < j < p, 
ij is a predecessor of ij+i. Such a sequence exists since rt G ^t+i. Let us prove by induction that 
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ij = rt — j: this is true for j = 0. Assume it is true until the integer j < p. We have either 
ij-^-i = ij — 1 or ij+i = ij + D — 1. But from the induction hypothesis, ij + D — 1 is an ancestor of 
ij+i, thus ij+i = ij — 1. This gives that s^_^_l is the largest column number i such that i < rt and 

^ ^ •5[i,t]- 

Now if we assume, by contradiction, that p > D — 1, then rt — D + 1 is not a predecessor of rj, 

which yields that rt has no predecessor, which is a contradiction. This gives the inequality of the 

first item. The second item is obvious, since s^^^ ha a unique predecessor which is s^j^i — D + l.D 

Lemma [T] induces a partition of fired columns between those which make a progress {i.e. in- 
creases the greatest fired column) and those which do not. This distinction is important in further 
development, so let us give progress firings a name. Let s^ = (sf,...,s^) be an avalanche, a 
column s^ is called a peak if and only if s^ > maxs^ ^_-^y 

Remark that two peaks p ^ q can be compared using chronological (<t) or spatial (<s') orders. 
Nevertheless, by definition of peaks we obviously have p <t q <^=^ p <s q. 

2.3. Pseudo local process 

The next Lemma explains precisely the way peaks appear, as soon as a D — 1 successive columns 
are fired. It follows an intuitive idea: a peak at time t + 1 is a column which only receives grains 
from the left part of the sand pile (within sm ^i)- Therefore, the amount it receives is at most 1 
and a peak must have an initial value of Z? — 1 units of height difference. Also, a non-peak column 
isn't fired when it receives 1 unit of height difference so it has to wait for its right neighbor to be 
fired, in a kind of chain reaction. 

Lemma 2. Let s'^ be the k^^ avalanche. Assume that there exists a column I such that p; / + Z? — 
2] C s , and a fired column i' G s such that i' > I + D — 1. Let I' be the lowest peak such that 
I' >l + D-l. 
There exists a time t such that [ if i is such that I' — D + l<i<l' then i ^ sf 



if i is such that I' < i then i ^ sf 



Moreover, let cr* denote the configuration obtained from iT{k — 1) via s^ ^, , then I' is the lowest 
integer such that V >l + D — 1 and at = D — 1. 



Informally, the lemma above claims that the space threshold I' induces a corresponding time 
threshold t: columns fired before time t are on the left of /', while columns fired after time t are 
on the right of /'. 

Proof. Let to be the time when s^^ = I', i.e. the first time such that s^^ > I + D — 1, and let j be 
the largest integer such that, for all / G [0; j], we have s^ _,_■, = s^^ — j'. Let us state t = to + j. 
We have j < D — 1. 

Let i, with i < I', such that i ^ s^ ^, . We claim that we have: i ^ s'^. To prove it, we prove by 

induction that for any t' > t, i ^ s^^ ^n. Assume that this is satisfied for a fixed t' . This means that 
all the transitions of sfj , ^^ ^'1 ^^^ done on columns larger than /'. Thus, aj = a\ and no transition 
is possible on i for cr* since s is leftmost (the only potential column to be fired is s^ — j — 1, but 
by assumption, either this column has been previously fired, or it cannot be fired by definition of 
j, according to Lemma [1]) . 
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By contraposition, it follows that for each column i G p; / + D — 1|, we have i G s^j^ ^, . A simple 

(reversed) induction shows that, for i G [/ + Z) — 1; /'J we have i G s^, since by hypothesis i + 1, 
and i + 1 — D both are in s . Thus, by contraposition of the claim above, for i G p + -D — 1; /'|, 
we have i G s^ ^i. This gives the the fact that for all i with /' — D + 1 < i < /', i G s^ ^ 

The fact that for all i with I' < i, i ^ s^^ n is trivial, by definition of to ai^d t. 

We have I' > I + D — 2 and (j\i = D — 1. assume that there exists /" < /' satisfying the same 
properties. Notice that for t^ < t' < t ^e have s^, > I' — D — 1. Thus the time ti such that 
•5*1 = ^" — -D + 1 is such that ti < to- That means that /" should have been fired before to, a 
contradiction. D 



Lemma [2] describes in a very simple way the behavior of avalanches. Thanks to it, the study of 
an avalanche can be turned into a pseudo linear execution, in which transitions are organized in a 
clear fashion: 



Theorem 1. Let s^ 



'1'' 



, Crpj 



he the k avalanche and {pi, . . . ,pq) be its sequence of peaks. 



Assume that there exists a column I, such that for each column i with i G p;/ + D — ![, i G s^. 
Then for any column p such that p > I + D — 1, 

p is a peak of s <;=^ 7r{k — l)p = D — 1 and 3j s.t. pj < p < Pj + D — 1. 

Furthermore, let j < q and t such that pj = s^, with pj > I + D — 1. Then 

T > t + Pj — Pj-i — 1 and for all t' s.t. t < t' < t + pj — Pj-i — 1, s^/ = s^/.i — 1. 

A graphical representation of this statement is given on figure [HI 



I 




X 


^>vAy VyV/ vvvyvAy 



I 

Figure 8: Illustration of Theorem [T] with D = 6; surrounded columns / to / + D — 2 are supposed 
to be fired; black column is the greatest peak strictly lower than 1 + D — 1; a column is grey if and 
only if its value is D — 1; following arrows depicts the avalanche 



Proof. The first part is a straight induction on Lemma [2l 

The second part follows from an induction summed up in the following fact: any column i such 

that 7r(/c — l)i<D — l must wait for its right neighbor z + 1 to be fired, and it should be fired when 

both i + 1 and i — D + 1 have been fired (besides, i — D + 1 has already been fired). Since any of 

such i is fired to reach a fixed point, we have for all peak pj > 1 + D — 1 that T > t + pj — pj^i — 1. 

D 



Theorem[T]easily allows us to compute the right part of the k avalanche (from column l+D—1 



only knowing 7r(A; — 1). The sequence of peaks is computed as follows. The first one is the lowest 
column i greater or equal to I + D — 1 such that -7r(/c — l)i = D — 1. Then, given a peak i, the next 
one is the lowest j such that 7r(/c — l)j = D — 1 and j — i < D — 1. If such a j does not exist, then 
there is no more peak and i is the largest fired column. 
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We can distinguish two movements within an avalanche: before a certain column I it has an 
unknown behavior, and from / to the end the behavior is pseudo local: when an index is fired 
ahead (on the right) then any 'hole' is filled before the progress can continue. Column / depends 
on parameter D and the number of grains N. We say that the k avalanche s is dense from I 
to m when m is the greatest fired column {\/ i > m, i ^ s'^) and any column between / and m 
included has been fired (V i G p; mj, i G s ). A consequence of Theorem[T]is that the avalanche s 
considered is dense starting at I, where I denotes the parameter in the statement of the Theorem. 
We will define the global density column C{D, N) as the minimal column such that for any avalanche 
s^, with k < N, s^ is dense from C{D,N). When parameters D and N are fixed, we sometimes 
simply denote C The formal definition of the global density column is: 

Definition 2. C'{D,k) is the minimal column such that the k^^ avalanche is dense starting at 
C'{D,k): 

C'{D, k) = min{/ e N | 3m G N such that V i G |/; mJ, i G s'' and \/ i> m,i ^ s^}. 

Then, the global density column C{D, N) is defined as: 

C(D, N) = max{ £'(£), k) \ k < N}. 

See figure [To] for an illustration of the global density column C{3,N). 

An important direct implication of Theorem [1] is that if there exists a column / such that for 
the k avalanche s , we have for all i G p; / + Z? — Ij, i G s , then for all j G p + D — 1; max s [, we 
have j — {D — l),j,j-\-l G s^. From the definition of the model, those three firings add respectively 
1,—D and D — 1 units of height difference to 7r{k)j. Since no other firing of the k avalanche 
affects TT{k)j, we have 7r(A:)j = 7r(/c — l)j. This behavior accounts for an important part of the fixed 
point, and this part keeps the same shape (see an example on figure [9]). But if the shape is the 
same, then it intuitively hints some similarity between successive avalanches. 



7r(1068) 5004155405 


2035545432154321 


;< > 1 1 1 1 1 1 


7r(1069) 0053054305 


20355454321054321 



£(6,1069) 

Figure 9: D = 6, 7r(1068) and 7r(1069) are represented as sequences of height difference. We 
pictured a part of §1069. g^ light grey square is a fired column and a dark grey square is a peak. 
By definition, s^o^g jg f^gj^ge starting from C{Q, 1069). We can notice that for every column j such 
that £(6,1069) + -D - 1 < j < 5^°^^, 7r(1068)j = 7r(1069)j, which hints that the next avalanche 
reaching that part of the configuration will behave in the same way. 



24. Application to KSPM(3): £(3, A^) in 0{\ogN) 

In this section we prove that in KSPM(3), the global density column £(3, A^) is in O(logA^). 
Considering the A^ first avalanches, from a logarithmic column in A" we can apply Theorem [1] and 
avalanches proceed pseudo locally, as described on figure El 
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Figure 10: D = 3, one avalanche is depicted on each hne. A grey square is a fired column. The 
black line illustrates the global density column C{3,N), and the dashed line depicts the maximal 
non-empty column. One can notice that the dashed line is 2 columns ahead of the rightmost fired 
column because it moves 1 grain to the new rightmost (maximal) non-empty column when applying 
the transition rule on that column). 

Proposition 2. [D = sjl £(3,iV) is in O(logiV). 

Proof. We prove that A^ is exponential in £(3, N). We use the fact that by definition of £(3, N) 
there exists an index N' < N such that £(3, A^) = £'(3, A'), which means that the A'*'^ avalanche 
is dense starting from £(3, A) but not dense starting from £(3, A) — 1, and then prove that A' is 
exponential in £'(3, A') (which proves the result since A > A'). For convenience and without loss 
of generality, let us consider that N = N' (to avoid using the symbol A^' and saying that A > A'). 
Let i = C{3,N), the trick is that for the particular value D = 3, the configuration 7r(A^ — 1) has 
unique and very regular values on columns 7r(A — l)o, vr(A — l)i, . . . , 7r(A — 1)^, 7r(A — l)i+i. 

We suppose £'(3, A^) > 0, which is true for A > 9. 

For convenience, let a = 7r(A — 1) and a = (oq, ai, . . . ) be its shot vector i.e., Oj is the number 



■^The proof given here is difTerent from the proof presented in [17 
of the present proof. 



We have more hope in the possible generahzation 
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/ 10 


We state A = 


1 




i^-1/2 3/2 



of times that column i has been fireqf|. According to the iteration rule we have the relation: 

(Tj = aj_2 — Soj + 2aj+i 
i.e., 

We denote by fj the column vector such that vf = (0,0, (Tj/2), 

and Ui the column vector such that uf = (aj_2,ai-i,aj) (with the convention that u is the line 
vector obtained by transposition of the column vector u) . The equality above can be algebraically 
written as 

Uj+i = Aui + Vi 

with the initial condition 

no = (A^-l,0,ao) 

Reduction of the dimension. We first get a simplification, exploiting the fact that the value 1 
is a double eigenvalue for A, as follows. Let (61,62,63) be the canonical basis of Z^, we consider 
the vectors: 

e'l = 61 + 62 + 63 

62 = 62 + 263 

We have Ae'^ = e'l and Ae2 = e'l + 62. On the other hand, we have 61 = e'l — e'2 + 63 and 
62 = 62 — 263 (which guarantees that (6'^, 62, 63) is a basis of Z^). 

Let P be the matrix of the projection on the vectorial 1-dimensional space Agg, generated by 
63, according to the direction of the plane Hg/ g/ , generated by e\ and e'2. We have: 




The plane Hg/ g/ is (globally) invariant by A, and, for each integer i, Ui — Pui is element of Hg/ 
Thus A{ui — Pui) is element of Hg/ g/ , which finally gives PA{ui — Pui) = 0, i.e., 

PAui = PAPui. 

So the equation lij+i = Aui + Vi gives Puj+i = PAu-i + Pvi, i.e. 

Pui+i = PAPui + Pvi 



•^For example, if we study £(3, 98) we can look at s^^ which is not dense starting from £(3, 98) — 1 and dense 
starting from £(3, 98), so £(3, 98) = £'(3, 98) (visible on figurefTU)). Then a = 7r(97) = (2, 0, 2, 0, 2, 1, 2, 2, 1, 0, 2, 1, 0") 
and a = (41, 14, 21, 12, 11, 7, 5, 3, 2, 1, 0"). Later in the proof we will have in this example that the sequence x — 
(138, -68, 34, -16, 8, -3, 2, 0, 1, 0, 0, 1, 0"). 
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We state Xi for the third component of Pui (both first and second components are null). The third 
component of Pvi is o"j/2. We have 

/ 
PA= 
\_i 1 _i 

Thus, the equality Puj+i = PAPui + Pvi yields to the equality: 

-Xi + (Jj 

with initial condition xq = N — 1 + uq. 

Relative expression of elements Xj. We can invert the previous equality to: 

Xi = -2Xi+l + CTj (1) 

We claim that the above inequality can be generalized to get: 

fc-i 
x,+i_fc = (-2)'^x,+i + ^(-2)'=-i-V,_, (2) 

r=0 

This is easy by induction. We have the initialization for k = 1, and if we assume the result for k, 
then we have 



Xi+i-{k+i) = -2xi+i_fc + ai_k = -2((-2)^Xi+i + ^(-2)^ ^ Vi_^) + a 



r=0 

fc-1 k 

x,.k = (-2)'=+^x,+i + ^(-2)'=-V,_, + ai.k = (-2)'+ia;,+i + ^(-2)'=-V,_, 

r=0 r=0 

which is the result. 

Specification in our case. We claim that /1'{3,N) is even i.e., /1'{3,N) = C{3,N) = 2j for 
some J E N and that a needs to have as a prefix the sequence 2,0,2,0,2,0,2,... until column 2j. 
That is, 2(02)-^ is a prefix of a. Indeed, the avalanche needs to propagate rightward until column 
2j (otherwise C'{3,N) < 2j) which is achieved thanks to values 2 on even columns, and all the 
odd columns 1, 3, . . . , 2j — 1 must not be fired (otherwise the hypothesis of Theorem[T]is completed 
and the avalanche becomes dense) which is completed by values becoming 2 when their right 
neighbor is fired. Any other value lets the avalanche become dense. 

Applying the general formula ([2]) with i = 2j and fc = 2j + 1, we obtain: 

2i 

X2, + i_(2, + l) = (-2)2^-+lx2, + l + j;(-2)2^-V2,-. 

r=0 

We have a2j-r = for r odd and cr2j_r = 2 for r even. Thus, stating 2s = 2j — r, we get 

xo = (-2)2^-+ix2,+i + 2 ^(-2)2- = -2(4)^X2,+i + 2^^^ (3) 

.=0 "^ 
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Conclusion. We now prove that xq is exponentially large in j. Using the fact that x, is an 
integer for all i {xi is the third component of Pui with Ui a vector of natural numbers and P an 
integer matrix), we consider 3 cases: 

1. |x2j+i| > 2. Then from equation ([T|), for any i <2j + 1 the sign of Xi is negative if and only 



4JH 



3 



if i is odd (signs are alternating and xq is positive), therefore from Q we have xq >2 

2. |x2j+i| = 0. Then equation ([3]) becomes xq = 2 ^ ^~^ . 

3. |a;2j+i| = 1. Since a"2J = 2, we have from equation ([1]) that either X2j = 4 if X2j+i = —1 
or X2j = if X2J+1 = 1. In the former possibility we get the sign alternation and apply the 
reasoning of case 1, and in the latter possibility we also have X2j-i = from equation ([1]) 
because o"2j-i = and apply the reasoning of case 2. 

In any case we get: 

4i _ 1 
N-l + ao = xo> 2—^ 

We obviously have oq < — g— since each firing of column moves 2 grains to the right. Thus, 
we get |(A^ — 1) > 2 7 from which we can easily deduce that: 

J < log4 {^N - -) 

hence, 

£(3,iV) = 2j<21og4(^iV_|) 

which gives the result. D 

For KSPM(3), after a short transient part of logarithmic length, the hypotheses of Theorem [T] 
are verified, and the study of avalanches can be turned into a pseudo linear process. Bounds on 
the maximal non-empty column size{D, N) of a fixed point with N grains shows that size{D, N) 
is in Q{\/N) when D is constant. Indeed any height difference 7T{N)i verifies < 7r(A^)j < D and 
it is easy to prove that there is no platearo of length greater than D columns, therefore size{D, N) 
is lower bounded by yjN/D and upper bounded by \J DN . As a consequence, the pseudo local 
process stands for the asymptotically complete behavior of avalanches. 

Unfortunately, the approach above does not hold for Z) > 3. The main reason is that, for 
D = 2> unfired columns induce a very particular and periodic prefix (2(02)-') on configurations. 
From D = A, the structure of such a possible prefix is more complex and we did not yet get a 
tractable characterization of it. 

3. Transduction 

In this section we study the temporal regularities — between successive avalanches — that can 
be derived from the density of avalanches. Therefore, the developments below hold starting from 
the global density column £(-D, N). 



*a plateau is a set of consecutive columns with the same height. 
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Firstly we study the similarities between two successive avalanches in subsection 13.11 We have 
seen in the previous section that when the k avalanche has been performed, a large part of the 
fixed point ir^k) is equal to 7r{k — 1), hinting a similarity between s^ and the next avalanche s^~^^. 
The similarity is stated in Lemma [3] by the mean of peaks: on a large part of the configuration, 
the peaks of s^ and s^^^ are exactly the same. 

Secondly, we use the equality between successive sequences of peaks (a kind of temporal stability 
of the peaks) and combine it with our description of the avalanche process (Theorem[T]) to construct 
a finite state word transducer in subsection 13.21 Given as an input a particular sequence of peaks: 
the greatest peaks between columns i{D — 1) and (i + 1){D — 1) for each of the A^ first avalanches, 
an iteration of the transducer outputs the greatest peaks between columns (z + 1){D — 1) and 
{i + 2){D — l) for each of the same A^ first avalanches. A second iteration of the transducer outputs 
the greatest peaks between columns {i + 2)(D — 1) and {i + 3){D — 1) for each of the A^ first 
avalanches again. The trace up to N on li = li{D — 1); (i + 1)(L' — 1) — 1] will be formally defined 
as the sequence of greatest peaks between columns i{D — 1) and (i + 1){D — 1) for each of the A^ 
first avalanches (a subtlety to be explained later is that a peak is counted once during its "lifetime" 
i.e., as long as it persists from an avalanche to the next one). Note that there may be more than 
one peak within columns of the interval Ij, but at most one peak of maximal value (the value of 
a peak is of course its index), which we name "greatest peak". We can therefore write that given 
the trace up to N on li, the transducer outputs the trace up to A^ on /j+i. 

Note that the transducer does not perform complex computation, because given A^ one can very 
easily compute the A^ first avalanches by starting from the empty configuration, adding a grain 
and performing the first avalanche, adding a grain and performing the second avalanche, and so on. 
This computation is done temporally: we compute s^, then s^, then s^, and so on. The interest of 
the transducer is that it embeds our knowledge on the regularity of avalanches and performs the 
computation spatially instead of temporally. Indeed, it computes the trace up to N on li, then on 
li+i, then on /j+2, and so on, each of the Ij representing a spatial "slice" of avalanches, as depicted 
on figure [TTl Now, the point is that if we know the trace up to A^ on Ij, we are able to compute 
7r(iV)fcforA:>(j + l)(Z?-l). 

The transducer depends on the parameter D of the model, as explained in the construction at 
the end of Lemma HI For D = 3 and A^ grains, we prove in subsection 13.31 (Corollary [T|) that for 
any trace up to N on li, there exists an index n in O(logA^) such that the trace up to A^ on Ij_(_„ 
is periodic. The proof uses only the transducer constructed for D = 3. Hence, we can take for i 
the smallest integer such that the transducer construction holds i.e., li is the first "slice" on the 
right of the global density column £(3, A^) (plus 1 "slice" for technical reason). Therefore i is in 
O(logA^) as proved in Proposition [21 and the periodic trace up to A'^ emerges on /j_|_„ with both i 
and n in 0(log A^). This study is concluded by Proposition O which states that when we compute 
7r{N)k ior k > {i + n + 1){D — 1) from a regular trace, it presents a regular wavy shape. We can 
conclude for D = 3 that the wave pattern emerges from a logarithmic column in the number of 
grain A^, thus the fixed points are asymptotically completely wavy. 

3.1. Successive avalanches 

Toward the construction of the finite state word transducer, we start with a formal description 
of the intuitive implication of Theorem [T] that two successive avalanches share a significant part of 
their peaks. Indeed, when the k avalanche is dense from / to m, for each column i G ll + D — 1; ?tt,[[, 
columns i — D + 1, i and i + 1 are fired within the k*^ avalanche. Therefore they both receive 
and give D units of height difference, and 7r(A:)j = 7r(A: — 1)| = 7r(A; — l)j (figured]). Moreover, 
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Figure 11: D=4. Avalanches reaching £(D, N) up to 500, one per hne. Since we exploit the density 
of every avalanche, an avalanche is meaningful for the study of the fixed point -7r(500) if it goes 
beyond the global density column £(4, 500) (and D — 1 = 3 columns further for technical reason). 
Non- meaningful avalanches are shadowed. The black line is the global density column £(4, A^), and 
the dashed line is the maximal non-empty column. A light grey square is a fired column, a dark 
grey square is a peak. The trace up to 500 on Is is 0120120210. From this trace, we can compute 
7r(500)fc for A;> 12 (12 is obtained as (3 + 1)(D - 1)). 

Each avalanche reaches a fixed point, so to every line A^ representing the avalanche s on the 
figure corresponds the fixed point tt{N), and to each position at ordinate A^ and abscissa k on the 
figure corresponds a difference of height iT{N)k- Following this convention, non empty columns of 
7r(500)fc for A; > 12 is hatched on the top right of the figure. 

A precise study of the transducer for D = 3 will lead to the prediction of a portion of fixed points 
(analogous to the one striped in this figure) which tends to 1 (subsection 13. 4p . 

For technical reasons, the trace up to A^ will be defined on slices on the left of column C{D, N) + 2{D — 1), therefore the trace 
up to 500 is formally not defined on Is. We nevertheless use this example for its meaningfulncss over size ratio. 
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TT{k)j = TT{k — l)j = for j > m + D — 1. An intuitive consequence is that two consecutive 
avalanches are similar. This intuition is formally stated in this section. 

Recah that s^ denotes the k^^ avalanche of KSPM(D). We define <^{D,N) = (c/)!, ..., (/>"), 
the subsequence of (s^, . . . , s ) reaching column C{D, N) -\- D — 1. Formally, s G ^{D, N) <;=^ 
C{D,N) + D - 1 G s''. ^{D,N) is called the sequence of long avalanches up to N of KSPM(i:)). 



An example is given on figure [12 
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Figure 12: Z? = 4. The 195 first avalanches are depicted on the left, only the long avalanches up to 



195 are depicted on the right — s € $(4, 195) 



£(4, 195) +D-le s''. The black line is the 



global density column £(4, N). A light grey square is a fired column, a dark grey square is a peak. 

We also define the corresponding sequence {fi^, fi^, , /i") of fixed points such that f/* = 7r(0) = 

0'^, and for each integer k, fi^ = 7r(m), where m satisfies 0'^ = s"^. 

The definition of long avalanche is motivated by Theorem [1] which states that the avalanche 
process is understandable in simple terms on the right of the global density column, therefore we 
concentrate on firings on the right of the global density column and thus consider only avalanches 
reaching that part of the configuration. 

Remark 1. Remark that in KSPM(L'), if s^ is a long avalanche up to A^, whose sequence of peaks 
is denoted by P (the largest peak being maxP ), from Theorem [1] we have (figure [9]): 

l)i foi C{D,N)+D- 1 <i <maxP''; 



7r(A;)i = TT{k - 


- 1 


'^{^)vaa.yiP^ = 


< 


■K{k.)i = 7r{k - 


-1 


ir{k)i = 7r{k - 


-1 



1) 



0; 



+ 1 for max P^ < i < max P^ + D 
for i > max P'^ + D — 1. 



1; 



This comes from a clear application of transition rules (for each considered column i we know 
which columns of the set {i — D + l,i,i-irl} are fired in s , so we can compute '7r(A;)j from 7r{k — l)). 
In other words. The knowledge of maxP'^ allows to compute n^k) from ir^k — 1) straightforwardly. 
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Lemma 3. In KSPM(D), let L he the global density column of N, and ^ = {cp^ , . . . , (p"^) its 
sequence of long avalanches up to N. Let k < N, and P'^ (resp. P^+^) be the sequence of peaks i 
of (j)^ (resp. (j)^~^^) such that i > L + 2{D — 1). We have: 

P'^'\{maxP'=} = P'^'+in|L + 2(D-l);maxP^[[. 

The lemma above can be seen as follows: |P^^| > \P \ — 1, and the |P | — 1 first elements of 
P^^ and P are equal. Informally, the peak sequence can increase in arbitrary manner, but can 
decrease only peak after peak. 

Proof. Let k, k! be two integers such that is the k avalanche, and </> "^"^ is the k! avalanche. 
For each column i such that z G \L + D — 1; maxP*^[[, we have i — D + l,i,i + l G 0'^ and therefore 
7r(K)j = ■k{k- l)i. 

By definition of long avalanches, any avalanche s between (jr and (p '^^ stops before L + D — 1, 
i.e. for all i>L + D — l,i^s. Combining it with previous remark, we have for all k" £ |k; k'J: 

for ahi G IL + D- l;L + 2(i:)-l)[[ , iriK"), > -k{k - l)i (4) 

for ahi G lL + 2(D-l);maxP*-'| , TT{K")i = tt{k - l)i (5) 

because columns within interval IL + D — 1;L + 2{D — !)[[ can gain height difference when a 
column within [L; L + D — If is fired. This is in particular true for k" = k' — 1. We now 
study the k'^^ avalanche (j)''^'^. From relation (JH) and since tt{k' — 1) is a fixed point, for all 
ielL + D-l;L + 2{D- !)[[, 7r(K - 1)./ = D-1 ^ tt{k' - 1), = D - 1. Let Q^ (resp. Q'^+i) be 
the sequence of peaks i of (j)'^ (resp. (j)^^^) such that «G {L + D — 1; L + 2{D — !)[, using Theorem 
[T]we therefore get 

Q'^ C Q'^+i (6) 

Let / = JL + 2{D - l);raayiP% From relation ^ 

for ah i G /, ■K{K-l)i = D-l ^=> Tr{K' - l)i = D - 1 (7) 

We now eventually prove the conclusion of the lemma. Let pi = min{i G P }, from Theorem [1] we 
equivalently have pi^ = min{z' G /|7r(«: — l)i = D — 1} (the existence of pj^ is a hypothesis of the 
lemma). Let p'-, = min{i' G P"^+i} = uim{i' G I\tt{k' — 1)^/ = D — 1} (the existence of p'-, is given 
by subset relation ([6]) ) , using relation ^ we have p'-, = pi. 
Other peaks within / are obviously equal from Theorem [1] and relation ([5]) with k" = k' — 1. D 

3. 2. Transducer 

We now exploit the similarity between successive avalanches. Intuitively, we will cut configu- 
rations into intervals Iq,Ii,I2, . . . of size D — 1, describe what happens in each interval via the 
notion of trace, and relate the trace on Ij to the trace on li^i with a finite state word transducer. 

The interval li is the column sequence {i{D — l),i{D — 1) + 1, ... , i(D — 1) + D — 2)). We call 
state of an interval /j of a fixed point vr its value (vrj(£)_i), 7rj(^_i')_j.X) • • • ) '^i{D-i)+D-2)- Hence, each 
interval state is an element of the set S = {0, 1, . . . , D — 1}^^^. An obvious, nevertheless important 
remark is that in 7r(0) the state of any interval is (0, 0, . . . , 0). It will enable us to compute a fixed 
point from a trace. For convenience, we also write 00 ... states and sequences when there is no 
ambiguity. 
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The formal definition of the trace up to N on an interval /, is more involved. Wordily, it 
is the ordered sequence of relative positions of the largest peaks p inside li for the N first long 
avalanches, where each peak is considered once while it stays from a long avalanche to the next 
one (this occurs when an avalanche performs firings strictly beyond p), see figure [TTl We formally 
define the trace in two steps: first the relative position of the largest peak inside /, — the type — 
then how we construct the trace by considering a peak only once. 

We fix i such that i{D — 1) > C{D, N) + 2(D — 1) and consider the interval /j. The maximal 
peak j of (/)'', such that j < {i + 1)(-D — 1), is denoted by p{i,k). The type a{i,k) of the long 
avalanche (j)'' on /j is defined as: 

• a{i, k) = p{i, k) mod [D — 1] if p(i, k) G /j; 

• a(i, A:) = e if p{i, k) ^ /j. 

Therefore, possible types are £,0,1,...,!? — 2. We formally consider types as words of length 
at most 1 over the alphabet T = {0, 1, . . . , D — 2}, and e as the empty word. So the set of types 
is TU {e}. Nevertheless, it is more natural to think about types as letters. 

Note that if a long avalanche (f)^ changes the state of /j , then from Remark [1] there necessarily 
exists a peak of (fr in the interval h-i- 

A trace is an element of T* . Intuitively, the trace up to N on /j is a subsequence of the whole 
sequence of types (a(i, A;))^^^, where we keep only once the value in each consecutive sequence of 
equal values. For example the sequence 01121002112220 should become 012102120. Toward this 
definition we introduce some more notions given an interval /j: type similarity between avalanches, 
subsequence of same type of avalanches, infiuent subsequence. 

We say that two long avalanches 0'^ and (j)^ are i-similar if they have the same type on /, i.e., 
if a{i,k) = a{i,k'). We can now divide the sequence $ of long avalanches up to N into maximal- 
length subsequences {(l)^,(/)^~^^, ....,(/)'' ) such that, for each integer k' S [[/c;A:"[[, (j)^ and cp'' ~^^ are 
z-similar. Such a subsequence is called an i- sub sequence. An i-subsequence is said of type a for i 
when the type of each avalanche of the subsequence is a. When a is not the empty word e, we say 
that the subsequence is i-influent. 

Definition 3. The trace up to N on li is defined as the sequence of type of the i-influent- 
subsequences of ^(D,N). 

Figure [13] illustrates the definition of a trace. 

Note that from Lemma El each (i + l)-influent-subsequence is contained in an i-influent- 
subsequence. Indeed, the peak sequence can be increased in an arbitrary manner but only the 
greatest peak can disappear from a long avalanche to the next one. Therefore if we consider an 
i-subsequence of type q, ^[k,k"] = {4>^-, ■ ■ ■ ,4''' ); ^[k,k"] ends only if the greatest peak of (j)'' is 
at relative position a in /,. Also, for the i-subsequence ^[k,k"] to start, a peak must appear on 
relative position a in /j during (j)^, and it enforces the previous avalanche <p^~^ to stop before this 
position (max(/) < i{D — 1) + q) because from Remark[T]the only changes induced by (j) on 
height differences which can create a peak are on the D — 1 columns following its largest peak. As 
a consequence 4)^~^ is of type e on Jj+i and so does not belong to an {i + l)-infiuent-subsequence. 
We can conclude that any (i + l)-influent-subsequence must start and end within an i-influent 
subsequence. 

The flnite state word transducer is constructed to compute the sequence of (i + l)-infiuent- 
subsequences within an i-influent-subsequence, hence relating the trace on /j to the trace on /j+i. 
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Figure 13: D = 4. Illustration of definitions around the trace up to 500 on 14. The left datas 
are extracted from figure [TTJ (T) Long avalanches up to 500, a dark grey square is a peak, a 
light grey square is fired column. £(4,500) = 6 and <1>(4, 500) = (0^, . . . ,4>^'^). @ Corresponding 
whole sequence of types on I^: (a(4, /c))^^^^. (3) 4-similar long avalanches are grouped into 4- 
subsequences. (3) The trace up to 500 on I4, is obtain with the 4-influent-subsequences: 0120120. 



Lemma 4. Let ^\k^k"] = {(j)'^,(j)^'^^, ....,4>^ ) be an i-subsequence of type a, with k" < N, and 
with li^i an interval whose columns are greater than C{D,N) + 3{D — 1). Given the state 
{ao,ai, ...,a£)-2) of li^i in the configuration fJ,^~^, and a, one can compute, with no need of more 
knowledge: 

• the state (aQ,a'^, ..., a^g) ''/-^j+i ^'^ ^^c configuration fi^ , 

• the sequence of types of the successive (i + l)-infiuent subsequences contained in ^\k,k"]- 

Proof. This is obvious when the type of the subsequence is e, since there is no change and the 
(i + l)-subsequence contained in {(p^ , ■ ■ ■ ,4>'' ) is also e. 

The computation is simple when there is no integer m £ [0; a] such that a^ = D — 1. In 
this case, the peak p{i, k) is the last peak of (j)^ , thus ^K^ ^^ = 0, which gives that p{i, k) is not a 
peak of (j)^^^ , thus the subsequence is reduced to a singleton which is not (i + l)-influent (second 



part of the result). For [D — l)(i + 1) < i < p{i,k) + D — 1, we have fi"j 



k-l 



^,- 



+ 1, and for 



p{i, k) + D-l < j < {D-l){i + 2), we have /i*- = ^^ . Thus, we have a[ 



, + 1 for m e [0; a] 



and a'^ = Om for m £}a; D — 2] (first part of the result). 

Otherwise, (p^ contains a peak in /j+i. Let q{i + l,k) denote the largest one. The column 
q{i + l,k) is the largest j such that /i'^ = D — I and j < {D — l)(i + 2). Thus q{i + l,k) 
mod D — 1 is the largest m such that Om = D — 1. In this case, (p'' starts an (i + l)-subsequence of 
type q{i + l, k). Consider the following long avalanches. From LemmaO while q{i + l, k) remains a 
peak of (j)^ , p{i, k) also remains a peak of (f)^ . From Remark[Tl while q{i + 1, k) is not the last peak 
of (f) , the state of /j+i remains invariant. So the first avalanche (fr that changes the state of Ij+i 
is the one whose last peak is q{i + 1, k) {k' < k" since ^[k^k"] is a complete i-subsequence). We have 
l^t(i+i k) ~ ^' which closes the {i + l)-subsequence of type q{i + l,k). We also have /i^ = ^^ + 1 for 



22 



q{i + l,k) < j < {D — l)(i + 2), and ^u^ = /i^ for p{i, k) < j < q{i + 1, k). This gives the state of 
/j+i for /i (as in the previous case, this can be rewritten to show that this state can be expressed 
only from a and (ao,ai, ...,0^-2)) and proves that p{i,k) =p{i,k' + 1). 

The argument above can be repeated as long as we have a column j of /j whose current value 
is D — 1. When there is no more such column, the peak p(i, k) is deleted (its value becomes 0) by 
the next long avalanche which is necessarily (p from the maximality of i-similar subsequences. D 

The algorithm below gives the exact computation underlying the proof Lemma [H From the 
state of an interval /j+i and an avalanche type on /j, / returns the greatest fired peak in /j+i, and 
g computes the new state of /j+i and appends the result of / to a sequence of types on interval 
/j_l_i. g recursively calls itself, anticipating the i-similarity of successive avalanches when maxP 
lies on the right of interval !{. 

Input: a non empty type a and an interval state A = (oq, • • • , aD-2)- 

Data structure: a sequence u of types. 

Functions: 



/:5xr^ru{e} 


g:SxTxT* ^SxT* 




/(Aa): = 

if {{m<a\am = D-l]^%) 

then 

max{m Cm = D — 1} 
else 

e 


g{A,a,u) := 
match f{A,a) with 

e -^ (ao + 1,... ,aa + l,aQ,+i,... 

\p -^ giiao, • • • , Op-i, 0, Op+i + 1, 


,0,0-2), u) 

. . ,aD-2 + l),a,u::p) 



Computation: (j4, a) 1— )• giA, a, e) 

The algorithm above allows to define a deterministic finite state transducer C (see for example 
[^) computes the trace up to N on /j-|_i given the trace up to N on li. It is 

Definition 4. The finite state word transducer C for D is a a 5-tuple ((5,S,r,/, 5) where: 

• the set of states Q is S; 

• the input and output alphabets (resp. S and T) are equal to A = T\ {e} = {0, . . . ,D — 2}; 

• the transition function 6 has type Q x S — )■ Q x F* and is defined by the algorithm above: 
6iA,a) = C omputationiA, a) ; 

• the initial state is (0,0, ... ,0), and we do not need to define a final state. 

The image of a word ti by C is denoted by tin). The formal definition of t will not be useful 
and is a bit involved, we nevertheless give it for the sake of formalness: 

Definition 5. In order to define t, we need to define some other functions. Let 63 (resp. Sf) 
denote the first (resp. second) projection of the result of the application of 5: 

if6iA,a) = iA',u), then ( 6siA,a) = A' 

\ 6riA,a) = u' 

We also define the generalized transition function of the transducer: 

5*iA,u) = 5°iA,u,e) with / 5°iA,au,u') = 5°i5siA,a),u,u' SriA.a)) 

6°iA,e,u') = iA,u') 
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Then t is defined for all word u £ A* by 



t{u) = <5f (00 . . . 0, n) 

where 6^ is the second projection of the result of the application of 5* . 

As presented in LemmalU if a is the type of an i-influent-subsequence ^r^ ^/i, and A the state of 
/j+i in jj!^^^, then 5{A, a) computes the state of Ij+i in p)' and the sequence of types of the {i + 1)- 
influent-subsequences within ^[k^m- It follows that if (3 is the type of the i-influent-subsequence 



$ 



[k'+l,k' 



1 , we already know the state of Ij+i in /i'^ and are able to compute the sequence of types 
of the {i + l)-influent-subsequences within ^[k'+i^k"]- 

We have already highlighted the fact that, in the fixed point 7r(0), any interval Ij is in the state 
(0, 0, . . . , 0). As a consequence, if u is the trace up to N on li, then t{u) is the trace up to N on 
/j+i. Furthermore, since the input and output alphabets of C are equal, we have that t{t{u)) is 
the trace up to N on Ij_|_2 and more generally that t"(u) is the trace up to A^ on /j+n- This has 
been illustrated on figure [6al where u is the trace up to N on C{D, N) + 3{D — 1). 

Note that, more generally, the last subsequence of a trace may not be completed. This detail 
will be discussed in subsection! 




Figure 14: Transducer for D = 3 - Edges are labelled x\u, where x G ^ is the type to the current 
interval (input) and u G ^* is the resulting sequence of types applied to the next interval (output). 
00 is the initial state. For example, t{abaaaaab) = abaab. Remark that, for n > 0, we have: 
t{{ab)'') = (a6)"-i. 

For the lowest interesting value, D = 3, the transducer C can easily be drawn. This diagram is 
given on figure [Hi For readability, we write a (resp. b) instead of (resp. 1) for the type alphabet, 
and we omit the drawing of states which are not connected with the initial one and are not useful 
for the computation of t{u), for any word u. 
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This transducer has three transient states, (00, 10 and 20) and four recurrent states (11, 12, 21 
and 22) organized in a cycle. A non trivial analysis of this transducer is given in subsection 13.41 
The result is stated on the following statement: 

Corollary 1. [D = 3] For any k there exists n in 0{logk) such that for all u of length k, t'^(u) is 
a prefix of (ab)^ . 

For D = 3, we therefore establish the emergence of a periodic trace consisting in an alternation 
of a and b. Note that for any N, the length of the trace up to A^ on any interval (where it is 
defined) is smaller or equal than N. Consequently, if we consider the trace up to N on the leftmost 
interval Ij such that i{D — 1) > C{D, N) + 3(1? — 1), Corollary [1] states that the trace up to N on 
/j_i_0(iogAf) is a prefix of (ab)^. Since D is fixed equal to 3 and recalling Proposition [21 we therefore 
have, for D = 3, that the trace up to N on Io{iogN) is periodic (figure [6b]) . The next subsection 
concentrates on the meaning of regular traces. 

3.3. From traces to waves 

Corollary [1] states that when the parameter of the model D is equal to 3, then a periodic 
trace emerges very quickly (on a logarithmic interval in the number N of grains) . The difficulties 
involved in the generalization of the emergence of regular traces to any parameter D are discussed 
in the proof of Corollary [TJ Simulations nevertheless let us believe that this behavior generalizes 
to any parameter, so we choose to present a generalized interpretation of regular traces. 

Proposition 3. In KSPM(D), let li he an interval whose columns are greater than C{D, N) + 
3{D — 1). Assume that the trace up to N on li is 

{0,...,D-2f{0,...,p), withx>0 andp< D-2. 

Lety be the length of the last i-influent- sub sequence of type p. We have y < x+1, and ^T{N)wi_^_lW£)_l\, 

equals 

{p,...,l){D-l,...,ir-yO{D-l,...,l)yO'- ify<x + l 
(p + l,...,l)(D-l,...,l)-0- ify = x + l. 

Proof. It is a straight induction on avalanches. We concentrate on the right part of fixed points: 
'''"(^)[(i+i)(D-i),oo[- Initially for A: = 0, it is equal to 0*^. The D — 1 first i-influent-subsequences lead 
to D - 1, D - 2, . . . , 1, 0"^. From (D - 1, . . . , 1)^0'^, the trace cyclically applies types 0,...,D-2. 
From Lemma[3l we can predict that each type will be repeated x + 1 times (there are x peaks to be 
consumed for each i-infiuent-subsequence) , thus each cycle consists of {D—l){x+l) long avalanches. 
Each cycle verifies the following invariant: the ((x + l)p + yY^ long avalanche, < p < D — 1 and 
< y < a; + 1, has type p and lead to 

p,p-l,...,l,{D-l,..., ir-yo{D - 1, ... , l)yo^ if y < x; 
p+l,p,...,l,iD-l,..., 1)^0'^ if y = x + 1. 

This invariant is proved at each step by a direct application of Theorem [TJ n 

Remark 2. Note that the application of Theorem [1] gives the trace up to N on /i+i, it is 
{0,. . . ,D — 2)^'~^(0, . . . ,p). As a consequence, (0, . . . , I? — 2)'^ is fixed point for t. 
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A regular trace for a parameter D is a prefix of(0,...,D — 2)^^ and 
a regular trace on an interval li implies a very regular and wavy shape 
starting from the next interval. The name "wave" is explained on figure 

m 




3.4- Analysis of the transducer for D = 3. 

In this subsection we provide an analysis of the transducer for D = 3, 
leading to a proof of Corollary [TJ Note that though we consider maximal 
length subsequences of long avalanches, input words for the transducer 
may contain arbitrary numbers of successive occurrences of a and b since 
we consider only i-influent subsequences. 

We need some notations. Let A and A' be states of S and w be a word 
of A*. Consider, in the transducer, the path which starts in A, whose 
sequence of successive edge (left) labels is given by u. We say that we have 
Au = A' if this path terminates in A' . A word u is an entry word if 00 u 
is a recurrent state and for each prefix u' of n, 00 u' is a transient state. 
Let tA '■ -A* ^ A* be the transduction function obtained by changing the 
initial state by A in the automaton. Hence too = t. We extensively use 

^21, so we state t2i = t' . A word u is basic for the state A ii \tA{u)\ > 2 and for each prefix u' of 
u, \tA{u')\ < 2. For each current state A, the set of basic words for A and their images by tA are 
given below (tables represent case disjunctions according the beginning of u) 



Figure 15: Regular 
traces imply reg- 
ular shape of the 
fixed point, visibly 
wavy. 
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Each word u (such that t{u) / e) admits a unique decomposition u = uoUi...Up such that uq is 
an entry word, for 1 < i < p, n, is a basic word for the state OOuoUi...Ui-i, and Up is a non-empty 
prefix of a basic word (for the state OOuoUi...Up^i). The word u also admits a decomposition 
u = u'iu'2---u' , such that for 1 < i < p, u[ is a basic word for the state 21u'qu'i...u[_^, and u' , is a 
non-empty prefix of a basic word (for the state 21u'qu'i...u' i_-^). 

A first result gives us a hint on the form of the sequence of types applied to successive intervals: 

Lemma 5. Let C be the language C = {abu, u £ A*} U {e, a}. 

• For each u £ A* , we have t'{u) G C 

• For each v £ C, we have t{v) £ L. 

• For each u £ A* , we have t^{u) £ C. 

Proof. We prove the three items successively, using previous ones as hypothesis. 
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• Let u £ A* such that u / e. Consider the second decomposition seen above: u = n'^Ug-.-u', . 
We obtain t'(n) = t'{u'i)tA{u2...u' ,), where A denotes a recurrent state. 

— For p' >2, t'{u[) is the image of a basic word for 21, thus t'(u'i) G {ab, aba}, which gives 
t'{u) G C. 

— For p' = 1, t'{u) = t'{u[) and t'{u'i) is the image of non empty prefix a basic word for 
21, thus t'{u[) is a prefix of aba, which gives t'{u) G C. 

• Let V £ C. If t; G {e, a}, then t{v) = e. Otherwise v can be written abu. Thus t{v) = 
t{ab)t'{u) = t'{u), and t'{v) G C from the first item. This proves: t{v) G £. 

• Let u £ A* such that n / e. We consider the first decomposition above: u = uoUi...Up. We 
obtain t{u) = tA{ui)tAui{u2---Up), where A denotes a recurrent state. 

— For p = 0, t{u) = e, thus t^{u) = e. 

— For p = 1, t(n) = tAl^i), and tyi(^i) is the image by Ia of a prefix of basic word for A, 
which gives that t{u) is a prefix of either aba or ba (since possible images of basic words 
are ab,ba, and aba). This gives that t'^{u) G {e, a}- 

— Ifp > 2, then tyi(^i) G {at>, ba, aba}. If tA(^i) G {ab,aba}, then t(ti) G £, thus t'^{u) G £, 
from the second item. If i^l^i) = ba, then we can state t{u) = bau' . Thus t^(n) = t'{u'). 
We have t'(w) G £ from the first item, thus t^{u) G £. 

D 

Definition 6 (Height). The height h of a finite word u £ A* is h{u) = \\u\a — \y\b\ where \u\x is 
the number of occurrences of the letter x in u. 

Lemma 6. For any finite word v £ C, we have: h(t{v)) < —j^ + 1. 

Proof. This is obvious if w G {e,a}. Thus, stating v = abu, it remains to prove that, for any 
finite word u G A*, we have: h{t'{u)) < — ^ + 1. 

Let us first consider the case when \u\a — \u\i, > 0. Assume that we remove a pattern of the 
form ab or ba from u. This does not change the value of h{u). Moreover, for each recurrent state 
A, tA{oh) and tA{ba) both are elements of {ab,ba} and Aab = Aba = A. This guarantees that 
pattern suppression does not change the value of h{t'(u)). 

Iterating this argument until there is no more pattern as above leads to the following fact: if 
we state u' = a^^'^\ then we have h{t'{u')) = h{t'{u)). 

The integer h{u) can be written as h{u) = Ai + r, with r G [[0;3]]. We have: t'{aaaa) = aba, 
and 21aaaa = 21. Thus t'{u') = {abay t' {a'' ) , which gives h{t'{u')) < h{{abay) + h{t{r)) <i + l. 
Thus h{t'{u)) < ^ + 1. 

The other case, when \u\a — \u\h < 0, is similar. By simplifications of factors ba and ab, we 
obtain that h{t'{u')) = h{t'{u)), for u' = b^^'^K The value h{u) can be written as h{u) = 4j + s, 
with s G [0;3]. We have: t'{bbbb) = abbab and 216666 = 21. Thus t'{u') = {abbaby t'{b'), which 
gives h{t'{u')) < h{{abbaby) + h{t{s)) = j. Thus h{t'{u)) < ^ + 1. D 

Corollary 1. Given a word u £ A* of length I, there exists an n{l) in 0{logl) such that t"('-'(n) 
is a prefix of (ab)^ . 
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Proof. We first prove it restricting ourselves on words of £ Given a finite word v on C, we 
define the maximal height g{v) = uia.x{\h{v')\ \v' prefix of v}. The previous lemma gives the result 
5'(i(i')) < 1 + 4^- We can now use a trick to get the expected result. We define g'{v) = g{v) — |, 
then: 

From lemma O t{v) is element of C. Thus we can iterate the inequality. By this way, we obtain, 
for each positive integer n: 



g'it'^iv)) < 



4?! 



Thus, for n > log4^{g'{v)) — log4(|), we have: g'{t"-{v)) < |, so g{t"'{v)) < 2 and, by integrity, 

This last inequality enforces that u admits a decomposition i"(t;) = wiW2---Wp such that, for 
i G |l;p], Wi e {ab,ba}, and Wp G {e,a,b}. Thus t"+^(n) = t{wi)t'{w2) ■ ■ .t'{wp). Thus, t"'+^(u) is 
a prefix of the infinite word (ab)'^ , since t'{ab) = t'{ba) = ab and t{ab) = t{ba) = e. 

Now, if we take a finite word u on A*, we have, from lemma O i^{u) G C On the other 
hand, \t^{u)\ < 4|n| and \t^{u)\ + | > g'{t^{u)), which gives g'{t^{u)) < A\u\ + |. Therefore, for 
n > log4(4|ii| + o) — log4(o), we obtain that t"+^(t^(n)) is a prefix of the infinite word {ab)^ . In 



3/ ^^&4V3; 

'"'' -^ 3^ -^^S4V3^ 



other words, for m > log4 (4|n| + o) — log4(o) + 3, t''^{u) is a prefix of the infinite word (ab)^. □ 



For D = 3, we can now prove the following result: 
Theorem 2. For D = 3 and all N, there exists a column n in 0{\ogN) such that 

vr(iV)[n,oo[ = (2,ir[0](2,irO- 

where * is the Kleene closure and [0] stands for at most one zero. 

Proof. Let A^ be given. From proposition [2l L{3,N) is in O(logA^), therefore there exists an 
index m in 0(log A^) such that traces up to N can be defined on Im, and the transducer for D = 3 
derived from Lemma [5] be applied. Corollary [1] tells that there exists an index / in O(logA^) such 
that the trace up to N on Im+i is a prefix of (ab)^ . Finally, Proposition [3] gives the result with 
n={m + l + l){D-l). D 

4. Conclusion 

This paper explored emergent regularities in the Kadanoff Sand Pile Model with parameter D 
in spite of its overall complex behavior. We first studied the sequence of avalanches triggered by 
the repeated addition of one grain on a stable configuration. Regularities in avalanches allow a 
precise description of their process by the mean of a distinguished set of columns, peaks. We then 
concentrate on a particular interval of columns of constant size, Ij, and keep track of any firing 
in this interval, while grain additions and avalanches are repeated. The obtained information, 
constituted by peaks indices is named a trace. The next step is the construction of a finite state 
word transducer which computes the trace on Ij_|_i from the trace on Ij. Applying again the 
transducer outputs the trace on Ij+2 , . . . The transducer is thus a finite tool which we can study 
in order to predict the behavior of avalanches, and hence the shape of fixed points. 

The main result of this paper is: 
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Theorem 2. For D = 3 and all N, there exists a column n in 0{logN) such that 

vr(iV)[n,oo[ = (2,ir[0](2,irO- 

where * is the Kleene closure and [0] stands for at most one zero. 

Computer aided simulations intimate that every resuh of our study generaUzes to any param- 
eter D. Nevertheless, some holes remain is this puzzle, summarized on figure [T6j We therefore 
conjecture: 

Conjecture 1. For a general parameter D and all N, there exists a column n in O(logA^) such 
that 

vr(iV)[n,oo[ = P - 1, D - 2, . . . , 2, 1)* [0](D - 1, Z) - 2, . . . , 2, 1)*0". 

Our approach gives a possible way for proving this conjecture. First, prove that avalanche den- 
sity is reached in 0(log N) columns. Second, understand the structure of the introduced transducer 
(and its iteration) for a general value D. 

N 



0{\ogN) O(logAf) 



e{^/N) 



■k{N) 




"> columns 



Results 


D = 3 


For all D 


®: C{D,N) \nO{\ogN) 


Proposition [2] 


X 


(2): Build a transducer 


Theorem [1] 
Lemma S] 


Theorem [1] 
Lemma [H 


(3): Regular traces 


Corollary d] 


X 


d): Wave pattern 


Proposition [3] 


Proposition [3] 



Figure 16: The precise study of the case D = 3 gives an asymptotically complete characterization 
of its fixed points according to the number of grains. Though experimentally confirmed, some 
parts of this study do not easily generalize to any parameter D of the model. We summarize in 
this picture our contribution, its results and conjectures. 

The left curve labelled denotes the global density column C{D,N) which has been proved in 
Proposition [5] to be in O(logA^) for the case D = 3. From C{D,N), we built a finite state word 
transducer to compute traces from interval to interval in Lemma HI using the avalanche process 
description of Theorem [TJ This transducer applies in the part (2), and Corollary [1] states that for 
D = 3 the trace up to N is periodic from a column in O(logA^) (label (3)). Finally, Proposition [3] 
establishes the relation between regular traces and wave patterns, on part @. 
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5. Further ^vork 

The work presented in this paper has been continued and improved. A different approach, 
extending the hnear algebra analysis introduced in subsection 12 . 41 for D = 3, is presented in [19!]. It 
uses a non-trivial change of basis of the linear system so that its behavior becomes understandable 
in simple terms, eventually leading to a proof of Conjecture [TJ 
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